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ABSTRACT 

Measurements of interstellar Na I AA5890, 96 absorption lines in 18 ultraluminous infrared galaxies 
(ULIGs) have been combined with published Na I data, in order to reassess the dependence of galactic 
outflow speeds on starburst luminosity and galactic mass. The Doppler shifts reveal outflows of relatively 
{N| ' cool gas in 15 of 18 ULIGs with an average outflow speed at the line center of 330 ± 100 km s _1 . 

, The relation between outflow speed and star formation rate (SFR), defined by the distribution's upper 

envelope over four orders of magnitude in SFR, demonstrates that winds from more luminous starbursts 
accelerate interstellar gas to higher speeds as roughly v oc SFR 35 . This result is surprising since, in 
the traditional model for starburst-driven winds, these relatively cool gas clouds are accelerated by the 
ram pressure of a hot, supernova-heated wind that exhibits weak (if any) X-ray temperature variation 
with increasing galactic mass. The lack of evidence for much hotter winds is partly a sensitivity issue; 
but the Na I velocities in ultraluminous starbursts actually are consistent with acceleration by the tepid 
' wind indicating a hotter component is unlikely to dominate the momentum flux. The Na I velocities in 

the dwarf starburst winds do not reach the terminal velocity of a hot wind at the measured temperature 
of kT ^0.73 keV - a result which could be interpreted simply as evidence that the hot superbubbles 
are too confined in dwarf starbursts to generate a free-flowing wind. A dynamically-motivated scenario, 
however, is that the dwarf starburst winds simply lack enough momentum to accelerate the clouds to the 
velocity of the hot wind. Among the subsample of starbursts with well-constrained dynamical masses, 
the terminal outflow velocities are found to always approach the galactic escape velocity. Motivated by a 
similar scaling relation for stellar winds, the galactic Eddington luminosity for dusty starbursts is shown 
to be within the range measured for ULIGs. If radiation pressure on dust grains, coupled to the cool 
wind, is indeed important for galactic wind dynamics, then feedback will be stronger in massive galaxies 
than previously thought helping shape the high-mass end of the galaxy luminosity function. Regardless 
of the nature of the acceleration mechanism in ULIGs, the mass flux of cool gas estimated from these 
data demonstrates that starburst-driven winds transport significant gas during the assembly stage of 
field elliptical galaxies - a factor which helps explain the rapid decline in SFR in these systems inferred 
. £^ | from elemental abundance ratios. 

Subject headings: galaxies: formation — galaxies: evolution — galaxies: fundamental parameters — 
ISM: kinematics and dynamics ISM: structure ISM: evolution 
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1. INTRODUCTION 

Galactic-scale gaseous outflows are ubiquitous in star- 
burst galaxies at all cosmic epochs (Heckman, Armus, 
& Miley 1990; Pettini etal. 2001; Shapley etal. 2003). 
Powered by supernova explosions, perhaps with additional 
help from AGN in very luminous galaxies, these galac- 
tic winds transport heavy elements, interstellar dust, and 
energy into galaxy halos and the intergalactic medium. 
The significance for galaxy evolution depends on details 
of the coupling between the available energy, both radia- 
tive and mechanical, and the interstellar gas. Observa- 
tions, which describe this interplay empirically, illuminate 
the dominant physical processes and provide quantitative 



constraints for galaxy evolution models. 

This paper directly addresses feedback at the high mass 
end of the galaxy luminosity function via measurements 
of outflow speeds in ultraluminous starbursts. These new 
observations of Na I and K I (ionization potential of 5.1 eV 
and 4.3 eV respectively) interstellar absorption reveal rel- 
atively cold gas clouds advected into the outflow at the 
shear interface between the hot wind and the quiescent 
gas disk. Material evaporated from these clouds could be 
the source of the mass-loading required to explain the large 
X-ray surface brightnesses of galactic winds (e.g. Suchkov 
et al. 1994; Suchkov et al. 1996; Strickland et al. 2002; Mar- 
tin, Kobulnicky, & Heckman 2002). Unlike emission lines, 
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possible by the generous financial support of the W.M. Keck Foundation. 
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absorption lines leave no ambiguity between infall and 
outflow and are sensitive to extended, low density gas. 
The results reported here complete a larger effort to de- 
scribe the low-ionization absorption kinematics over four 
orders of magnitude in star formation rate (Heckman et al. 
2000, luminous infrared galaxies (LIGs); Schwartz & Mar- 
tin 2004, dwarf starbursts). The observations have little 
overlap with the Rupke et al. (2002) sample of ULIGs due 
to different selection criteria (see §2). 

Feedback is expected to have a large impact on dwarf 
galaxies owing to their shallow gravitational potential. It 
can explain the absence of young galaxies with very low 
mass, V ro t^, 10 — 20 km s" 1 , after reionization (Barkana 
& Loeb 1999), the slope of the galaxy luminosity func- 
tion (Dekel & Silk 1986; Dekel & Woo 2003; Benson etal. 
2003), and the mass-metallicity relation among galaxies 
(Larson 1974). Ablation of the gas in dwarf galaxies by 
winds from nearby galaxies may help flatten the faint-end 
of the luminosity function (Scannapicco, Thacker, & Davis 

2001) . 

Observations strongly support the idea that a larger 
fraction of the metals are removed from dwarf galaxies 
than massive galaxies. Metallicity measurements from X- 
ray spectra indicate the hot wind carries the heavy ele- 
ments produced by the starburst population (Martin et al. 

2002) . Temperature estimates for the hot gas, i.e. kT x , 
vary little with galactic mass indicating a critical galaxy 
mass, V c <~ 130 km s _1 , below which most of the hot wind 
escapes (Martin 1999). Retention of a larger fraction of 
the metals in more massive galaxies would produce an in- 
creasing effective yield with galaxy mass up to this scale 
where the stellar yield is recovered (Larson 1974; Dekel 
& Silk 1986; Vader 1987). Although the X-ray data have 
some significant shortcomings - including a small number 
of observed objects and a temperature biased towards the 
densest regions of the hot plasma, the mass-metallicity re- 
lation has been measured locally for a complete sample 
of galaxies. The metallicities increase steeply with mass 
over the range from 10 8 5 - 10 10 5 MqHjq but flatten at 
10 10 - 5 M Q ft,f 2 (Trcmonti etal. 2004). The mass scale of 
the flattening appears to be consistent with the inferred 
mass scale of metal retention (Garnett 2002). Despite the 
efficiency of metal ejection, it is not clear whether the to- 
tal mass loss rates which are of the order of the SFR are 
sufficient to explain the measured faint-end slope of the 
galaxy luminosity function, which is flatter than that of 
the halo mass function (Somerville & Primack 1999). 

This paper finds an increase in the velocities of cold 
clouds with galactic mass suggesting feedback may play a 
more prominent role in massive galaxies than previously 
thought, an idea which could be important for several out- 
standing problems. First, wind velocities are quite im- 
portant in determining how much of the IGM is polluted 
(Aguirre 2001). Fast winds from quasars may well deposit 
more energy into the IGM than starburst winds (Scan- 
nepieco & Oh 2004), but they seem unlikely to dominate 
the metal transport. Second, the galaxy luminosity func- 
tion cuts off more steeply than the theoretical halo mass 
function at bright luminosities, and a feedback mechanism 



physically distinct from that in dwarf galaxies appears to 
be required to explain both ends of the luminosity func- 
tion. The mass scale separating these forms of feedback 
shows up not only in the mass-metallicity relation and 
X-ray temperatures but also characterizes the galaxy mass 
scale ~3x 10 10 M Q where (1) surface-brightness stops ris- 
ing with luminosity and mass and (2) less massive galax- 
ies have stellar populations weighted towards younger ages 
(Blanton etal. 2003; Kauffmann etal. 2003a, b). Any feed- 
back processes invoked to explain these relations must do 
so without disrupting the Tully-Fisher and Faber- Jackson 
relations between velocity, i.e. the depth of the grav- 
itational potential, and stellar mass, V oc M X J A (e.g. 
Bcrnardi etal. 2003), which show no obvious change at 
this scale (Zwaan etal. 1995; Sprayberry etal. 1995; Dale 
etal. 1999). 

Sample selection and the observational strategy are de- 
scribed in §2. The kinematics and column densities of 
the nuclear ULIG spectra are presented in §3. Section 4 
develops scaling relations between outflow properties and 
galaxy parameters, including the dynamical state of the 
merger. The results are shown to challenge the standard 
model of supernova-driven winds in §5, and the role of ra- 
diative acceleration is discussed. §6 summarizes the main 
results. The second paper in this series (Paper II) de- 
scribes the surprisingly large spatial extent of Na I ab- 
sorption, compares kinematics of Na I and Ha lines, and 
discusses the outflowing mass flux in detail. 

2. OBSERVATIONS 

The spectra presented here extend Na I studies of dwarf 
starbursts (Schwartz & Martin 2004) and luminous in- 
frared galaxies (Heckman etal. 2000) to ultraluminous 
galaxies (ULIGs). Targets were chosen from the IRAS 
(60 micron) 2 Jy sample (Strauss etal. 1992) to have 
bolomctric luminosities greater than 5 x 10 11 L Q , a 60^tm 
bump F' 2 (60/im) > i^(12/mi)x.Fi / (25/zrn) ) and declination 
8 > —35° . Spectral observations of 41 of the 64 galaxies 
satisfying these criteria were obtained. The selection was 
random unlike the Rupke et al. (2002) sample which pre- 
selected galaxies with outflows from the 1 Jy survey. The 
rather large redshift range, z = 0.042 to z = 0.16, reflects 
the paucity of ultraluminous galaxies locally. 

Echellete spectra were obtained 2000 September 19-20 
and 2001 March 26-29 using ESI on Keck II. The nucleus, 
identified in R and K band images (Murphy etal. 1996), 
was placed on the slit. The position angle was usually 
chosen to include the second nucleus. If the object could 
not be observed when the parallactic angle swept through 
the position angle defined by the two nuclei, then an an- 
gle closer to parallactic was chosen. The median seeing 
was 0'.'8 FWHM (full width at half-maximum intensity), 
which resolves the galactic disks. The 1" slit provided a 
resolution of ~ 70 km s _1 . 

The data were reduced using the echelle package in 
IRAF. 4 Fixed pattern noise was removed from the CCD 
frames using bias frames and internal quartz illumination 
frames. Fitted arc lamp exposures (Cu+Ar+Xe) provided 
a dispersion solution accurate to ~ 0.15 A. The nuclear 



4 IRAF, the Image Reduction and Analysis Facility, is written and supported by the IRAF programming group at the National Optical Astron- 
omy Observatories (NOAO). NOAO is operated by the Association of Universities for Research in Astronomy (AURA), Inc. under cooperative 
agreement with the National Science Foundation. 
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spectra were traced and extracted, and a sky spectrum 
extracted from the same slit was interactively scaled and 
subtracted. The continuum was fitted and normalized to 
unity in each order. Velocities quoted in the paper are 
converted to the rest frame and given relative to the Local 
Standard of Rest. 

The primary uncertainty in the outflow velocity of the 
cool gas is the systemic velocity of the galaxy. This paper 
focuses on the CO-sample which contains 15 ULIGs with 
CO velocities from Solomon et al. (1997) and 3 ULIGs 
with CO velocities from Dr. Aaron Evans (pvt. comm.). 
The stellar velocities derived from the Mg I lines are fully 
consistent with the CO velocities although the uncertain- 
ties were typically larger. The CO emission was chosen 
mainly for its insensitivity to extinction and the follow- 
ing dynamical consideration. In galaxy - galaxy mergers, 
the interstellar gas sinks toward the dynamical center of 
the merged galaxies as the orbital angular momentum is 
shed at large radii. The central gas surface density grows 
to more than 10 times the maximum in the Milky Way, 
and the ISM is almost entirely molecular (Sanders etal. 
1988a). The measured CO velocity therefore provides a 
better description of the systemic velocity than cither HI 
21-cm emission, which comes from a larger region, or emis- 
sion lines from HII regions. For the CO-selected sample 
described in this paper, the typical accuracy of the CO ve- 
locities is ~ 20 km s" 1 (Sanders etal. 1988b). In the CO- 
sample, the minimum luminosity is 6.5 x 10 11 Lq, and the 
implied star formation rates are several hundred M© yr _1 
as shown in Table 1. Four of the galaxies are members 
of the well-studied Bright Galaxy Sample (Sanders et al. 
1988b; Sanders etal. 2003). 

3. RESULTS 

Figure 1 shows the Na I absorption lines in the nuclear 
spectra. Three features are immediately obvious. First, 
the doublet lines which are separated by 300 km s _1 are 
blended together in all the spectra, so the line widths are 
much larger than those observed in dwarf starbursts and 
LIRGs. Second, the total Na I equivalent- widths, see Ta- 
ble 2 are large - ranging from 1 to 10 A. And, third, 
the lines are not black at the line center. In this section, 
we determine the contribution of interstellar absorption to 
these line profiles, describe the kinematics of the absorbing 
clouds, and estimate the Na I column density. 

3.1. Distinguishing Interstellar and Stellar Na I 

Interstellar absorption lines in galactic spectra are 
formed against a continuum which contains photospheric 
features. Since the alkali metals are mostly neutral in the 
atmospheres of late-type stars, such stars present promi- 
nent Na I absorption lines. The strength of excited pho- 
tospheric features, which have no interstellar counterpart, 
such as the Mg I A5167.32, A5172.68, A5183.60 triplet 
constrain the contribution of stellar Na I absorption how- 
ever. In ESI spectra of AFGK dwarf and giant stars, taken 
in the same manner as the galaxy data, the strength of the 
NaD doublet relative to the Mg I triplet is well-described 
by the empirical relation EWnuD — ^/SEWMgi- 

In galaxy spectra, the Mg I lines are blended with 
other photospheric lines, particularly an Fe I line, because 
they are broadened by 200-300 km s" 1 FWHM. Figure 2 



compares the total Na I equivalent width to the mea- 
sured equivalent width of the entire Mg I+Fe I complex, 
rest-frame bandpass from A5162 to A5189, in the nuclear 
ULIG spectra. Excited photospheric features are strongest 
in IRAS00153+5454, IRAS00188-0856, IRAS17208-0014, 
IRAS20087-0308, and IRAS23365+3604; but these galax- 
ies also present Na I absorption much stronger than that 
predicted by the stellar locus. Although the Na I strength 
for IRAS00262+4251 exceeds the predicted stellar equiv- 
alent width by only one standard deviation of the mea- 
surement errors, the significant Dopplcr shift of the lines 
suggest an interstellar component is present. In Figure2, 
the upper limit for IRAS 19458+0944 is not very restric- 
tive, but the upper limit for the other null detection, 
IRAS16487+5447, allows only a stellar component. 

Table 2 lists the total NaD equivalent width (in the rest- 
frame A5874 to A5899 bandpass) and the fraction inferred 
to be stellar from the scaling between stellar Na I and 
Mg I. The stellar contribution to the NaD equivalent width 
is typically ~ 10% and will be ignored in the remainder of 
this paper. 

3.2. Interstellar Gas Kinematics 

Since the stellar Na I absorption is negligible, the line 
profiles describe the kinematics of interstellar gas clouds 
and, if a resolution element is optically thin, the Na I 
column density at a particular velocity. Dopplcr shifted 
absorption components are a robust signature of outflow 
(for a blueshift) since the absorbing gas must lie on the 
near side of the galaxy. One obtains no direct information 
about the location of the absorbing gas along the sight- 
line, however. As discussed in §2, the systemic velocities 
for the CO-subsample are well determined. The system- 
atic analysis of the Na I absorption presented here draws 
the robust conclusion that nearly all ULIGs drive fast out- 
flows. For any particular ULIG, the interpretation may be 
complicated by absorption from tidal debris (e.g. Norman 
etal. 1996; Hibbard, Vacca, & Yun 2000). Tidal debris, 
in contrast to starburst-driven outflows, subtend a small 
solid angle (L. Hernquist, pvt comm), so they are unlikely 
to be the dominant source of the blue-shifted absorption 
lines. 

3.2.1. Doppler Velocities 

The spectral region around Na I AA5890, 96 is plotted 
in Figure 1. Comparison of the systemic velocity, marked 
by a doublet, to the line profiles reveals that much of the 
absorbing gas is Doppler shifted to wavelengths shorter 
than the that of the rest-frame doublet. To describe this 
dynamic component, a second doublet was fitted while the 
velocity of the systemic component was held fixed. Al- 
though the best fits, which are optically thick, are sum- 
marized in Table 3, another set of models with the ratio of 
the line strengths held equal to the ratio of the oscillator 
strengths (two-to-one) is provided in Table 4. Comparison 
of the velocities in these two tables illustrates the maxi- 
mum systematic uncertainties in the measurements. 

Although the fully resolved line profiles might be more 
complex than these simple Gaussian fits assume, the fit- 
ted velocities of the dynamic absorption component, vb 
in Tables 3 and 4, provide a fair estimate of the mass- 
weighted outflow velocity. The average fitted outflow ve- 
locity is 330 ± 100 km s _1 . The uncertainty in ^ is 
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the same neutral gas. The doublet spacing is 1327 km s _1 
so these lines are not blended. Figure 3 shows that the K I 
absorption profile closely follow the shape of the Na I pro- 
file in IRAS10565+2448 and IRAS20087-0308. Although 
the K I absorption profile for IRAS15245+1019 reveals an 
absorption line blueward of the Na I terminal velocity, this 
candidate high velocity component is not detected in the 
K I A7699 line (not pictured) and therefore is not the K I 
A7696 line. The K I lines were not detected in the nuclear 
spectrum for IRAS17208-0014, and residuals from subtrac- 
tion of strong OH sky lines prevent an accurate measure- 
ment of the K I terminal velocity in IRAS18368+3549. 
Hence, when the K I doublet was detected, it presented 
the same terminal velocity as Na I. 

In summary, among the 18 galaxies in the CO- 
subsample, 15 present neutral gas outflows. The 
galaxies IRAS16487+5447 and IRAS19458+0944 show 
no (or weak) interstellar Na I absorption. Only the 
IRAS08030+5243 line profile shows a rcdshiftcd compo- 
nent. 7 This outflow fraction of at least 80% in ULIGs, 
which is higher than the outflow fraction in luminous in- 
frared galaxies and (Heckman et al. 2000) dwarf starbursts 
(Schwartz & Martin 2004). The most luminous LIGs also 
presented a higher outflow fraction than lower luminosity 
LIGs, but the significance of the trend was unclear because 
the two subsamples were selected differently. The new 
ULIG data establish a trend of increasing outflow fraction 
with luminosity. Indeed, all ULIGs may have outflows, 
and this hypothesis is explored further in §4.3. 

3.2.2. Line Widths 

In Table 3, the average line width of the dynamic compo- 
nent is < Av B >= 320± 120 km s" 1 FWHM, much larger 
than the thermal velocity dispersion of warm neutral gas 
in galaxies. For comparison, an average sightline through 
the Milky Way disk detects 6 clouds per kpc in Na I - 
each with velocity dispersion ~ 8 km s _1 (Spitzer 1968). 
Those LIGs with interstellar-dominated Na I absorption 
also present very broad lines (Heckman et al. 2000). While 
turbulent motion, particularly at the interfaces between 
the hot wind and the disk (e.g. Figure 11 in Heckman 
etal. 2000), is expected to broaden the lines, numerical 
simulations of bipolar outflows appear to have difficulty 
producing the large observed velocity range (D. Strick- 
land, pvt comm). 

To produce a width of 300-400 km s _1 , a sightline 
may need to intersect multiple turbulent sheets or shells, 
each with distinct bulk motion. Higher resolution spec- 
tra (6 km s _1 in Figure 1 of Schwartz & Martin 2004) 
support this claim. Along a sightline into the center of 
M82, the broad Na I doublet, 236 km s" 1 FWHM, is re- 
solved into 5 distinct velocity components, each of width 
25 - 79 km s _1 . And, although the Na I profile remains 
smooth in NGC 1614, the more optically thin K I lines do 

5 When D 3 ys is varied from unity to its optically thin limit, D aya = 1.98, the systemic component absorbs more flux in the center of the total 
line profile; and, as a consequence, the velocity of the dynamic component shifts to larger outflow velocities. 

6 In the IRAS00262+4251 fit, the fit statistic is not acceptable for the unsaturated model, so the -293 km s _1 of the saturated model best 
describes the outflow. In IRAS 10494+4424, the saturated model describes the two minima better than the unsaturated model. 

7 The absorbing gas in IRAS08030+5243 appears to be falling toward the CO-defined systemic velocity at 250 km s _1 . Since no redshifted 
flows are detected in the other 17 galaxies, and an outflow from one galaxy can be detected in absorption against the continuum from another 
galaxy (see Paper II), one wonders if the orientation of the sightline toward IRAS08030+5243 is unusual or special. The object 10" southeast 
of IRAS08030+5243, however, has been classified as a foreground star based on its surface brightness profile, so there is no obvious companion 
galaxy (Murphy etal. 1996). 



only a few km s _1 when the profiles show multiple local 
minima. For example, superposition of the shorter wave- 
length line from the systemic component and the longer 
wavelength member of the dynamic component creates 
three local minima in IRAS08030+5243 and IRAS17208- 
0014. Two minima and a strong blueshifted line wing 
require two blueshifted components in IRAS10565+2448 
and IRAS18368+3549, and three components were fitted 
with one component held at the systemic velocity. Fitting 
smooth, blended line profiles, on the other hand, does not 
yield a unique solution, and these systems have uncertain- 
ties up to ±70 km s _1 in in Table 3. Selection of a 
low doublet ratio for the systemic component systemati- 
cally lowers the fitted velocity of the dynamic component. 5 
Comparison of vb in the optically thin limit, Table 4, 
and the optically thick limit, Table 3, shows the magni- 
tude of the shift is generally negligible. Two exceptions 
are IRAS00262+4251, where the shift is 100 km s" 1 , and 
IRAS 10494+4424, where the discrepancy is 110 km s _1 . 
In both cases, the model in Table 4, which has the satu- 
rated systemic component, is favored. 6 Systematic uncer- 
tainties in vb caused by the choice of the doublet ratio are 
therefore insignificant compared to the magnitude of the 
outflow velocities. 

The maximum outflow velocities were measured where 
the Na I line profile intersects the continuum. These veloc- 
ities presumably reflect the terminal cloud velocity signif- 
icant distances from the starburst. These measurements 
are not influenced by the details of the line fitting, so they 
are tabulated with the basic measurements in Table 2. 
The mean terminal velocity is -750 km s _1 , which is a 
few times larger than the galactic rotation speeds. This 
coincidence is particularly interesting since the escape ve- 
locity from the galactic halo is also a few times the rotation 
speed. The physical meaning of this correlation is explored 
further in §5. 

The He I 5876 emission - present in all spectra except 
IRAS00188-0856, IRAS03158+4227, IRAS18368+3549, 
and IRAS20087-0308 - could cause the true terminal ve- 
locity to be systematically underestimated. At velocities 
greater than 729 km s _1 , absorption from the 5889 Na 
line is observed at the same wavelength as systemic He I 
5876 emission. Three factors indicate this coincidence does 
not bias the measured terminal velocity. First, the appar- 
ent strength of the He I 5876 emission line relative to the 
upper limits on the He I 4471 line (another triplet with the 
same lower state as He I 5876) and He I 6678 (the n=3 
to n=2 singlet transition) indicates little, if any, of the 
He I 5876 emission is cancelled by very high velocity Na D 
absorption. Second, the Na I absorption extends further 
along the slit than the He I emission in many cases, and 
Paper II shows similar terminal velocities are measured 
all along the slit. Finally, the K I AA7665, 99 doublet, al- 
though a weaker line than Na I, would be expected to trace 
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break up into multiple components. An important test of 
dynamical models is to determine, presumably using nu- 
merical simulations, whether a more realistic spatial dis- 
tribution of star formation can generate multiple outflow 
shells and a broad enough velocity distribution. 

The fitted systemic Na I component probes a mixture 
of stellar absorption and interstellar absorption at the sys- 
temic velocity. The formal uncertainties for the widths of 
the systemic component are significantly larger than those 
for the dynamic component because the systemic compo- 
nent is generally weak. The average width of the systemic 
component Av sys is 400 ±140 km s _1 FWHM corresponds 
to a line-of-sight velocity dispersion, a = 170 ± 60 km s _1 
. These widths are substantially larger than those mea- 
sured for stellar-dominated LIGs (Hcckman ct al. 2000) 
but lie at the low end of the range measured in ULIG 
spectra from stellar lines <r* = 140 — 290 km s _1 (Gen- 
zel etal. 2001). The Na I sample presented in this pa- 
per has three galaxies in common with the Genzel et al. 
sample. In IRAS17208-0014, the stellar velocity disper- 
sion, <r* = 229 ± 15 km s _1 and v* = 110 km s _1 , is 
subtantially larger than that of the systemic Na I com- 



ponent, a b 



FWHM/2.35 = 108 km s" 



How- 



ever, in IRAS23365+3604 the systemic Na Icomponent, 
cr S j /s (Na I) w 259 km s _1 is much larger than the stel- 
lar velocity dispersion, cr* = 145 ± 15 km s _1 and 
= 15 km s _1 . The widths are similar in IRAS20087- 
0308 - a* = 219 ± 14 km s" 1 , v* = 50 km s" 1 , and 
(Na I) = 203 km s" 1 ; but this may be largely coinci- 



the Na I lines are not guaranteed to be optically thin. The 
columns in Table 4 are good lower limits. The fitted dou- 
blet ratios, D = £W(A5890)/£W(A5896), in Table 3 il- 
lustrate how much larger the true Na I column might be. 9 
The fit statistic was improved for every galaxy when the 
doublet ratio was allowed to decrease from the optically 
thin limit, D = 1.98. Figure 4 shows that the thick models 
do a better job of fitting the structure in the line profiles. 
Optically thin fits cannot be ruled out for the smooth, fea- 
tureless profiles in IRAS00262+4251 and IRAS20087-0308, 
but there is no obvious reason to expect these systems to 
be optically thin when the better constrained fits all indi- 
cate optically thick outflows toward the nucleus. 

3.3.2. Cloud Covering Factor 

Optically thick clouds that cover the continuum source 
would leave absorption lines completely black at the line 
centers. None of the lines in Figure 1 are black at line cen- 
ter. The conclusion that the Na I lines are saturated in 
many, if not all, of these galaxies directly implies that the 
absorbing clouds do not completely cover the continuum 
source. 

In the limit of large optical depth, the cloud cover- 
ing factor Cf is simply w 1 — /5890, where /5890 is the 
residual intensity in the A5890 line. The covering factors 
listed in Table 2, which use the /5890 value fitted in Fig- 
ure 1, indicate that 25% of the visible continuum source 
is covered by the clouds on average. The covering fac- 
tors show a relatively large range, however, from 10% in 
IRAS16090-0139 and IRAS00262+4251 up to nearly 70% 
in IRAS15245+1019. 

When optically thick gas only partially covers the con- 
tinuum source, the measured equivalent width will be in- 
sensitive to the Na I column density and will be primarily 
determined by the product of the covering factor and the 
line-of-sight velocity dispersion. Figure 5 shows that the 
FWHM of the dynamic component is not correlated with 
the total Na I equivalent width. It does not correlate with 
the EW of the dynamic component either. The bottom 
panel of Figure 5 shows that the lowest equivalent width 
systems all present low covering factors, and that large 
covering factors are only found in galaxies presenting very 
high Na I equivalent width. A large covering factor of cool, 
outflowing gas may partially reflect the viewing angle, but 
the analysis of §4.4 indicates temporal variations are also 
important. 

3.3.3. Conclusions about Column Densities 

Seven of the ULIGs show extremely strong NaD equiv- 
alent width (EW > 7 A); they are IRAS00188-0856, 
IRAS08030+5243, IRAS10565+2448, IRAS15245+1019, 
IRAS17208-0014, IRAS18368+3549, and IRAS20087- 
0308. The optically thin models generally give larger Na I 
columns in the galaxies with larger Na I equivalent widths. 
The association is not exactly one-to-one however. For ex- 
ample is IRAS15245+1019 is 5th in total EW but 3rd in 

8 For gas and dust mixed in roughly the Galactic ratio, sightlines with iV(Na I) > 1.7 X 10 22 cm -2 would have visual extinction Ay > 10 
completely obscuring the continuum source. 

9 The models adopt an optically thick systemic absorption component. The strength of the A5896 line in the systemic component is limited 
by the overall line profile, so the effect of increasing the doublet ratio from 1.0 to 1.98 is to strengthen the A5890 line. Since the Dopplcr shift 
of the dynamic component superimposes the A5896 line of the dynamic component on the systemic A5890 line, the fitted A5896 line of the 
dynamic component becomes weaker and the doublet ratio of the dynamic component is larger (i.e. lower column). 



u sys \ 

dental considering that the uncertainty in Av sys exceeds 
100 km s _1 . It is not surprising that the fitted Na I sys- 
temic components do not appear to be accurate measures 
of the stellar velocity dispersion - their strength suggests a 
blend of stellar and interstellar absorption at the systemic 
velocity. 

3.3. Na I Columns 

Measurements of the Na I column density in the dy- 
namic component provide a basis for estimating the mass 
flux of outflowing material per unit area. These measure- 
ments differ from absorption measurements using quasars 
in that the continuum is extended (slit width subtends 
^1 — 2 kpc of galaxy) and extremely dusty sightlines are 
excluded from the average. 8 In principle, measurements of 
the K I lines are more sensitive to column density than the 
Na I lines, which are shown to suffer from saturation; but 
subtraction of bright sky lines near K I generally leaves 
large uncertainties about the K I line strength. 

3.3.1. Fitted Models 

The inferred Na I column is only directly proportional to 
the measured Na I equivalent width if the absorbing mate- 
rial is optically thin. Assuming optically thin clouds, the 
derived Na I columns range from ~ 10 12 cm~ 2 to a few 
times 10 13 cm" 2 , see Table 4. Since the optical depth of 
the stronger doublet member, Na I A5890, reaches unity at 
a column of A(Na I) w 3.6 x 10 13 cm" 2 (v/100 kms _1 ) 2 , 
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total column because the particularly large Doppler shifts 
leave very little absorbed flux near the systemic veloc- 
ity. Hence, the fraction of the total equivalent width pro- 
vided by the systemic component is particularly low in 
this galaxy. The relative equivalent widths of the dynamic 
and systemic components can be ambiguous when the pro- 
file is smooth. Among the 7 largest equivalent width sys- 
tems, this happens in IRAS20087-0308 and IRAS00188- 
0856 where fits with more absorption at v sys are clearly 
allowed. 

Some of the profiles absolutely require very large fit- 
ted columns. Focusing on the systems with the largest 
equivalent width again, one finds that the systemic 
component has to be very weak in IRAS 10565-2448, 
IRASf8368+3549, and 1RAS15245+1019. Similarly, the 
systemic component is strong but maxed out in the 
fits to IRAS17208-0014 and IRAS08030+5243. In sum- 
mary, freeing the doublet ratio raises the Na I column to 
A(Na I) <~ 10 14 cm~ 2 to f 15 cm~ 2 in galaxies for which 
the optically thin model indicated large columns. The 
correction, however, has little impact on the columns in 
weaker systems, iV(Na I) x 10 12 cm~ 2 since it is not linear 
and therefore increases the range of measured Na 1 column 
densities. 

4. DISCUSSION 

The absorption measurements presented here trace gas 
clouds with a significant neutral component. As discussed 
above, the absorbing material is thought to be disk gas 
entrained in a hot supernovae heated outflow. The X-ray 
observations are biased by conditions in the highest den- 
sity regions of hot component and do not necessarily re- 
veal the hottest (lowest density) regions of the hot outflow 
(Strickland & Stevens 2000), but they do show that most 
of the shock heated mass can escape from dwarf galaxies 
with velocities less than 130 km s _1 (Martin 1999; Mar- 
tin 2004). More of the shock heated gas remains bound in 
larger galaxies, and this differential feedback contributes 
to the flattening of the slope of the faint-end of the lu- 
minosity function. No direct measurements of the kine- 
matics are currently possible in the hot phase. Since the 
clouds absorbing in Na I are presumably embedded in this 
medium, it is possible to learn about the wind fluid itself 
via the motions of the clouds. 

4.1. Outflow Mass 

The inferred H column associated with a sightlinc de- 
pends on the chemical abundance, dust depletion, and ion- 
ization state of the gas. Assuming solar abundances in the 
ULIGs, the H column is 

-- 4 - 9oxio20cm - 2 (i^) (¥ 

( AW \ 
\N(NaI) J 1 ' 

Heckman et al. (2000) argue, based on parameters of 
Galactic clouds, that the depletion of Na onto dust 



Observations of the absorption-line profiles at multiple 
positions across the galactic disks of ULIGs, presented in 
Paper II, indicate that some outflows cover more than 
10 kpc of the host galaxy. Taking a fiducial radius of 
~ 5 kpc for purposes of illustration here, the mass of 
clouds, M ~ NnAfimp or 



M r ~ 1.08 x 10 8 M m 



N Na 



10 14 cm 
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are applied. The kinetic energy carried by the cold 
clouds, l/2Mv 2 , is 



E ~ 1.71 x 10 56 ergs 



N. 
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The enormity of these quantities is perhaps most appreci- 
ated when the mass loss rate, M c = Sir 2 pv = Qr/rnifiNHV 
or 



M c = 141 M Q yr- 1 



N H 



R 



4.9 x 10 20 cm- 2 J \5 kpc 



400 km 



is compared to the SFR. For example, applying these ex- 
pressions to IRAS10565+2448 indicates a mass flux in cool 
gas of ~ 300 M yr -1 , which is similar to the SFR of 
190 M Q yr -1 . Paper II presents spatially resolved line 
profiles and more detailed estimates of the mass loss rates, 
but this example illustrates the potential importance of the 
cold outflows. If they can be sustained for even 10 Myr, 
several billion solar masses of interstellar gas is transported 
out of the disk thereby starving the starburst for fuel. 

4.2. Variation of Outflow Velocity with Galactic Mass 

One of the principal motivations for this study and that 
of Schwartz & Martin (2003) was to investigate how out- 
flow velocities scale with galactic mass. Surprisingly, Heck- 
man et al. (2000) found no correlation between outflow 
velocity and galactic rotation speed, or starburst lumi- 
nosity. Simple extrapolation of this result would imply 
dwarf galaxies completely dominate intergalactic enrich- 
ment. High-resolution Na I spectra of seven dwarf star- 
burst galaxies, however, revealed outflows in three sytems; 
but the velocities, —25 km s" 1 to -35 km s -1 , were much 
lower than those in the LIGs (Schwartz & Martin 2004). 
Provided their luminosities arc dominated by the starburst 
rather than a quasar, the ULIGs have star formation rates 
grains, d^ a = (N(Na) gas + A r (A?'a) ( z MSt )/A r (A''a) sas , is from 190 - 750 M Q yr -1 , which (by definition) is higher 



likely a factor ~ 10, and the ionization correction, e.g. 
N(Na)/N(NaI), could be as much as a factor of 3 to 10. 

The nuclear sightlines sample a small region of each 
galactic halo, so extrapolating these column densities to 
masses is very sensitive to the chosen area. 



than that for LIGs which drop down to about 10 M© yr " . 
Outflow velocities, all measured in Na I absorption, can 
now be compared over four orders of magnitude in SFR. 

In Figure 6, these velocities exhibit scatter at a given 
SFR, but the upper envelope of the velocity distribution 



7 



clearly rises with SFR. The galaxies plotted in Figure 6 
are all selected for their starburst activity. They repre- 
sent the highest areal star formation rate at a given lu- 
minosity and, presumably, the strongest outflows. Ve- 
locity measurements in normal galaxies would fill in the 
low velocity regions of Figure 6. The observed velocity 
spread may reflect mainly projection effects, as argued in 
§4.3. In this linear least squares fit to these data, 

log(V) = (0.35 ± 0.06) log{SFR) + 1.56 ± 0.13 where the 
outflow velocity V is in km s _1 and the SFR has units 
M Q yr -1 , should describe the slope of the upper envelope 
in Figure 6. 

The dynamical masses of the galaxies are not as 
well constrained as the star formation rates but pro- 
vide a more fundamental quantity for comparing the out- 
flow velocities. The three dwarf galaxies have rotation 
speeds of 35 km s" 1 (NGC 1569; Stil & Israel 2002 ), 
40 km s- 1 (NGC 4214; Walter et al. 2001), and 87 km s" 1 
(NGC 4449; Hunter, vanWoerden, & Gallagher 1999) and 
are clearly less massive than the LIGs and ULIGs. Ul- 
traluminous activity, on the other hand, requires a large 
amount of gas, M gas ~ 10 10 M , and a major merger 
(Borne etal. 2000). The expectation of massive progeni- 
tor galaxies is confirmed by measured rotation curves and 
stellar velocity dispersions for three ULIGs in this sample 
(Genzel et al. 2001). In the LIG luminosity range, however, 
the SFR is a poor surrogate for dynamical mass. The LIG 
subsample for which Hcckman et al. had rotation curves 
has a mean rotation speed of 199 ± 45 km s" 1 , which is 
less than that of the Galaxy indicating some LIGs are not 
massive enough to become ULIGs. Some fraction of LIG 
progenitors, however, are likely major mergers that were 
ultraluminous near perigalactic passage but faded during 
the long journey through apogalacticon. This idea is sup- 
ported by (1) the significant overlap in CO masses be- 
tween LIGs and ULIGs (Gao & Solomon 1999), (2) the 
revived starburst activity suggested by the bimodel dis- 
tribution of ULIG merger ages (Murphy etal. 2001), and 
(3) the wider mean separation of LIGs relative to ULIGs. 
Unfortunately, most of the LIGs with measured rotation 
speeds are highly inclined disks for which no outflow was 
detected by Heckman etal. (2000). Consequently, rota- 
tion curves are only available for two of the LIGs shown in 
Figure 6. Considering the paucity of data and poorly de- 
termined inclination corrections for the outflow velocities, 
these results should be taken as preliminary Nonethe- 
less, the combined data do indicate a steep rise in out- 
flow velocity with rotation speed, illustrated in Figure 8. 
This result is potentially of fundamental significance, so it 
will be important to obtain more rotation curves. For the 
galaxies without rotation curves, rough circular velocities 
estimated from line-widths yield few velocities above than 
the preliminary relation, so Figure 8 may actually provide 
a reasonable description of the upper envelope. 

Absorption lines are ideal probes of the kinematics of 
extended low density gas, so it is not surprising that 
outflows were one of the first features identified in spec- 
tra of high-redshift galaxies, whose optical spectra in- 
clude the plethora of strong absorption lines in the rest- 
frame ultraviolet bandpass. The choice of systemic ve- 
locity is critical when comparing outflow speeds. In par- 
ticular, Lya emission is often seen redshifted relative to 



stellar features due to resonance scattering off the reced- 
ing side of the flow and therefore makes a poor zero- 
point. Among Lyman-Break-selected galaxies at z ~ 3, 
the galaxy MS1512-cB58 has the highest quality spectra 
due its large lensing amplification. The measured SFR, 
40 M Q yr _1 , and outflow velocity, 255 km s _1 (Pcttini 
et al. 2002) place it right on the average relation for nearby 
LIGs. The high-redshift objects with bolometric luminosi- 
ties in the ULIG range are selected in the sub-mm (Ivison 
etal. 2000). Rest- frame ultraviolet spectra show median 
offsets between the Ha velocity and low-ionization ultra- 
violet resonance absorption lines of ~ 650 km s _1 at lu- 
minosities of a few hundred M yr (S. Chapman, pvt 
comm), which place them in the same part of Figure 6 
as the local ULIGs. Outflows from high-redshift galax- 
ies therefore appear to obey the same scaling relations as 
winds from local starbursts. 

Neutral outflows are not only faster in more massive 
galaxies, but their terminal velocity also appears to in- 
crease almost linearly with the galactic rotation speed. 
The bottom panel of Figure 8 shows the terminal veloc- 
ities are always 2 to 3 times the rotation speed. This 
normalization is particularly interesting since the escape 
velocity is \[2 (minimum) to 3.5 (isothermal halo extend- 
ing to 100 kpc) times larger than the circular velocity. 
The estimated escape velocities are roughly 100 km s _1 , 
400 km s _1 , and 900 km s _1 in dwarf starbursts, LIGs, 
and ULIGs respectively The proximity of the upper enve- 
lope in Figure 6 to these values provides important infor- 
mation about the dynamics of the outflows and is discussed 
further in §5. 

4.3. Velocity Variations from Sightline Orientation 

X-ray and optical imaging of nearby starburst galax- 
ies indicate the outflow axis is aligned perpendicular to 
the disk plane. Identical bipolar outflows viewed at ran- 
dom angles will yield a large range of measured Na I 
velocities, so it is interesting to examine whether the 
fastest outflows might be found in systems with galactic 
disks oriented nearly face-on to our sightline. This hy- 
pothesis is supported by the paucity of outflows in the 
edge-on sample of LIGs (Heckman etal. 2000). Incli- 
nation is only well defined for a few of the ULIGs due 
to both the lack of high-resolution imaging and the on- 
going merger. NICMOS images of IRAS17208-0014 and 
IRAS10565+2448 show these disks are oriented close to 
face on (Scoville etal. 2000), and the outflow velocities in 
both are among the largest measured at comparable SFR. 
If large outflow velocities are purely an inclination effect, 
then ULIGs like IRAS23365+3604 and IRAS18368+3549 
must also be observed nearly face-on, which is not in- 
consistent with their appearance in ground-based images. 
Large axial ratios indicate two ULIGs are likely viewed 
edge-on. Of these IRAS03521+0028 has very weak Na I, 
and IRAS11506+1331 has little net outflowing material in 
Na I. Inclination therefore appears to have something to 
do with observed velocity spread at a given SFR. 

To examine inclination effects quantitatively, consider a 
simple model where the flow is perpendicular to the disk 
plane. The average polar viewing angle is i = 57° , and 
the average projected velocity is 0.5v max - It follows that, 
at any particular SFR, the true outflow velocity is twice as 
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large as the average projected outflow velocity measured. 
Multiplying the fitted vb vs SFR relation by 2.0 yields 
the dashed line sketched in Figure 6. This estimate of 
the de-projected outflow velocity describes the upper en- 
velope adequately for all but a few, particularly luminous, 
objects. The distribution of observed velocities about the 
mean is shown in Figure 7, where the data were divided 
by the mean velocity at the appropriate SFR and then 
binned according to their deviation from the mean. The 
naive, planar flow model predicts a constant number of 
galaxies with increasing v/v avg up to v max = 2.0u a „ g . The 
data show more galaxies with outflow speeds near the av- 
erage than predicted, but this may be consistent with a 
more realistic representation of the wind geometry where 
the velocity field has a radial component. Projection ef- 
fects, however, fail to explain the high velocity tail in the 
distribution, and this shortcoming is further discussed in 
§5- 

4.4. Merger Sequence 

In addition to inclination, the other property that 
clearly varies among the ULIGs is the dynamical age of 
the merger, which is constrained by morphology (Murphy 
etal. 2001). Visual and near infrared imaging of ULIGs 
drawn from the 2 Jy sample show galaxy pairs, tidal tails, 
and double nuclei (Murphy et al. 1996). Inspection of these 
images allows some ULIGs to be identified as first passage, 
late passage, and fully merged. The ESI spectra resolve 
a number of nuclei, previously classified as single nuclei, 
into double nuclei. The outflow velocities do not correlate 
with dynamical age, but the strength of the wind does 
show some interesting trends that may provide some in- 
sight about the role of the winds in the ULIG-QSO tran- 
sition. 

Highly eccentric orbits are required to create an in- 
tense, short-lived disturbance which centrally concentrates 
the interstellar gas, so the galaxies pass rapidly on their 
first encounter. Short tidal features are one signature of 
first passage because it takes time for tidal features to 
grow. Disk-like morphology also indicates not enough time 
has elapsed for the stellar orbits to reflect the transferred 
orbital energy Of the 18 ULIGs discussed in this pa- 
per, five are plausibly on their first passage. Short tidal 
features are seen in IRAS10565+2448, IRAS11506+1331, 
IRAS19458+0944, and IRAS17208-0014, and the debris in 
IRAS00153+5454 is confined to the region around the dis- 
torted disks. These are all double nuclei systems. Of the 
five strongest Na I lines (i.e. equivalent width > 7 A), two 
are found among these 5 galaxies. The other three strong 
lined systems are caught during the second pass. In Fig- 
ure 9, these two are pictured along with their line profiles 
to represent the first-pass mergers. 

The duty cycle of the ultraluminous phase must be 
< 20% because the estimated dynamical timescales for the 
mergers, 5 x 10 8 yr to 15 x 10 s yr (Murphy etal. 1996), 
exceed the gas consumption timescales, ~ 10 8 yr. Most 
of the time between passages is spent near apogalacticon 
where the galaxies are unlikely to be ultraluminous. None 
of the systems with strong winds are found in a widely 
separated phase. The lack of tails and largely symmet- 
ric structure in IRAS00188-0856, IRAS16487+5447, and 
IRAS 19297-0406 point toward a stage just before the sec- 



ond encounter. The relatively large separations 20.2 kpc, 
5.3 kpc, and 4.4 kpc are consistent with this picture. Their 
Na I absorption, Figure 9, is quite weak, comparable to 
Na I line strengths in LIGs. These observations indicate 
the wind strength decreases when the galaxies are far apart 
between encounters. Three systems in the ULIG subsam- 
ple can be identified with the second (or nth) passage by 
their double nuclei, large tidal structures, and/or more 
spheroid-like than disk-like structure. The identification as 
2nd-passage systems is most secure for IRAS15245+1019, 
IRAS20087-0308, and IRAS18368+3549, which have par- 
ticularly long tidal tails. All three of these ULIGS present 
very large Na I equivalent widths. 

The nuclei were not resolved in 7 systems, and the 
interstellar Na I lines are generally weaker in these sin- 
gle nuclei systems. The most recently merged sys- 
tems still present large, diffuse tidal structures as seen 
in IRAS00262+4251, IRAS03521+0028, IRAS16090-0139, 
and IRAS23365+3604. The line strengths are moderately 
strong falling between those of systems on the nth-passage 
and the most advanced mergers in this sample. The other 
two systems, IRAS03158+4227 and IRAS08030+5243, 
present quite symmetric and much fainter tidal debris; and 
the interstellar Na I absorption is weak. These results sug- 
gesting the wind is dying down in the older systems. 

The combined small size of this ULIG sample and the 
large variations in the accuracy of the dynamical ages esti- 
mates makes a robust statistical analysis unfeasible. How- 
ever, this preliminary exploration strongly suggests that 
the strongest winds are found in systems near perigalac- 
tic passage. Furthermore, the winds appear to weaken 
when (1) the galaxies are far apart and (2) the merger is 
complete. The starburst probably dies down due to grow- 
ing scarcity of gas. Further investigation may determine 
whether the gas is simply consumed, removed by the star- 
burst, or removed by a nascent AGN-driven outflow. 

5. ACCELERATION OF THE COOL WIND 

Consistency between the X-ray and Na I measurements 
is an important test of the standard dynamical model for 
starburst winds, which assumes a thermally-driven hot 
wind accelerates cooler, entrained gas (Chevalier & Clegg 
1985; Heckman, Armus, & Miley 1990). Numerical simu- 
lations based on this idea reveal that large vortices develop 
at the shear layer between the hot, bipolar wind and the 
gaseous galactic disk (e.g. see Figure 11 in Heckman etal. 
2000). The Na I absorption is thought to be a good tracer 
of the cool gas entrained along these interfaces. The mo- 
tion of the clouds is therefore expected to reflect the prop- 
erties of the hot wind fluid. The X-ray surface brightness is 
biased towards the densest part of the wind, however; and 
observations may not detect the hottest component of an 
outflow. Nonetheless, two factors motivate a simple model 
where the wind temperature does not vary much with star- 
burst luminosity or mass. First, the X-Ray temperature 
of the mass-carrying component of the hot wind does not 
show a significant dependence on galactic mass (Martin 
1999; Heckman etal. 2000). Second, theoretically, the hot 
wind temperature is simply the amount of thermalized su- 
pernova energy per unit mass of entrained material. The 
energy and entrained mass may both scale linearly with 
starburst luminosity producing similar temperatures in all 
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will be accelerated to a maximum velocity determined by 
the equation of motion. Provided the cloud lifetime ex- 
ceeds the outflow timescale, the ram pressure of the hot 
flow accelerates a cloud of mass M r and 
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hot wind streaming 
momentum injection 
which is about 5 x 
starburst (Leithcrcr 



where v w is the velocity of the 
out through solid angle f2. The 
rate is p w = rh w v w = Up w r 2 vl,, 
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etal. 1999). 11 The absorbing clouds must enter the hot 
wind at relatively low velocities because the absorption 
line profiles extend smoothly to the systemic velocity. Re- 
placing the time derivative in Equation 5 with vdv/dr, ig- 
noring the gravitational deceleration, and integrating out- 
ward from the sound speed of the cool clouds, i.e. ~ 6 
km/s, yields terminal cloud velocities that increase with 
SFR as illustrated in Figure 10. Up until the point where 
the clouds are moving with the hot wind, these solutions 
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starburst winds. Another model, especially applicable to 
the ULIGs, is that radiation pressure on dust grains in the 
clouds accelerates the cool wind. Similar ideas regarding 
this radiative feedback have been advanced for limiting the 
mass of OB associations (Scoville etal. 2001). This idea is 
particularly interesting since the associated Eddington-like 
luminosity for galaxies would limit the maximum bright- 
ness of a galaxy in a given gravitational potential (Scoville 
2003; Murray, Quataert, & Thompson 2004). In what 
follows I give some initial interpretation of the data pre- 
sented in this paper. A more extensive discussion is being 
developed in a series of papers with Eliot Quataert, Todd 
Thompson, & Norm Murray. 

5.1. Supernova-Driven Winds 

The rise in Na I outflow velocity with galactic mass, 
Figures 6 and 8, appears to be in direct conflict with 
the relatively uniform X-ray temperatures measured over 
a similar mass range (Martin 1999; Heckman etal. 2000; 
Martin 2004). Gas at T x ~ 8.9 x 10 6 K (0.76 keV) 
in a hot bubble allowed to stream through a rup- 
tured supcrshcll accelerates to a terminal velocity of 
VZc s w 500(^/0.76 keV) 1 / 2 km s" 1 in the intcrgalactic 
medium. This velocity is remarkably similar to the speed 
of the Na I clouds in the ultraluminous starbursts. The 
mystery is thus reduced to understanding why the cool 
gas is not accelerated to the wind velocity in less luminous 
starbursts. 

One can imagine several reasons why the clouds might 
not reach the wind velocity. It would be interesting, for 
example, to examine whether the gaseous halos of dwarf 
galaxies can prevent the formation of a free-flowing wind 
(cf. Silich & Tcnorio-Tagle 2001; Mac Low & Ferrara 
1999). The X-ray temperatures indicate that winds in 
massive galaxies and dwarf galaxies start out with similar 
specific thermal energy, but the observations do not have 
enough spectral resolution to directly measure the kine- 
matics of the hot fluid. 10 Alternatively, the clouds might 
be photoionized or disrupted via hydrodynamic instabil- 
ities before they reach the velocity of the hot wind, but 
this solution is not appealing, however, since it requires a 
complicated scaling of cloud lifetime with starburst lumi- 
nosity. A simpler hypothesis is that the hot wind in lower 
luminosity starbursts does not carry enough momentum 
to accelerate clouds to the hot wind velocity. If correct, 
then the cloud velocities should tend toward a constant 
luminosity in more luminous starbursts. 

The upper envelope of the velocity - SFR relation in 
Figure 6 may flatten beyond a SFR of > 10 M Q yr -1 . 
The data points in this region are sparse, but most would 
agree that a break in the slope is not ruled out. Since 
X-ray observations suggest the mass flux in the hot wind 
is similar to the SFR, or cquivalcntly, that the mass load- 
ing is about a factor of ten larger than the mass ejected by 
supernovae and stellar winds (Martin 1999; Martin etal. 
2002), the rate of momentum injection into the hot wind 
is relatively well constrained. A given column of cool gas 

10 An even hotter component may be present in massive galaxies, but current limits on its emission measure suggest it carries a tiny fraction of 
the mass. Indeed, as long as more powerful winds entrain more interstellar gas, it is hard to avoid similar mass-weighted X-ray temperatures 
among starbursts of all luminosities. 

11 The Starburst 99 population synthesis code was used to estimate the rate at which supernovae and stellar winds supply momentum, i.e. 
p* = m« je ji>», e j = y / 2m» ! ejii»- A Salpeter stellar mass function from 0.1 Mq to 100 Mq was adopted to be consistent with Figure 6, and the 
star formation rate was assumed to be constant for 100 Myr. At an age of 10 Myr, p w would be about a factor of two lower. 



for clouds of density p c launched from radius R (e.g. 
Murray etal. 2004, Appendix A). Substituting N c = 
p c p~ 1 mp 1 R c , the fiducial wind speed of 500 km s _1 , 
rh w — AM*, and an arbitrary launch radius of 200 
the maximum cloud velocity can be written 
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where A <~ 1 and p was taken to be 1.4. This approx- 
imation indicates the critical SFR required to accelerate 
a cool gas column of Nh ~ 10 20 , 10 21 , and , 10 22 cm~ 2 
to 500 km s" 1 is 0.24 M Q yr" 1 , 2.4 M© yr" 1 , and 
24 M© yr -1 , respectively. Taking proper account of the 
relative cloud - wind velocity in Figure 10 shows somewhat 
higher SFR's, M, ~ 0.7, 7, 70 M Q yr" 1 , are required to ac- 
tually flow with the wind. Inspection of Figure 10 indicates 
that any combination of launch radius and column density 
satisfying (^A/10 22 c m - 2 )(i? /200 pc) ~ 1 is a good de- 
scription of the upper envelope of the outflow velocities 
for all but a few extremely luminous ULIGs. The implied 
rate of momentum injection, p v 

p w = 8.0 x 10 34 dyne {n/ir)N 22 R^ 2Q0 v 2 o 

which, for any solid angle 0, is easily supplied by a SFR of 
70 M Q yr -1 since p\ = 3.5 x 10 35 dyne. Supernova-driven 
winds can, it turns out, describe the observed increase in 
Na I velocities from SFR 0.1 to roue hly 100 Mq yr" 1 and 
still be consistent with X-ray temperatures that are inde- 
pendent of galaxy mass. 
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It should be emphasized that the terminal cloud ve- 
locity is sensitive to the launch radius because most of 
the cloud acceleration happens at small radii where the 
ram pressure is largest. If the size of the starburst region, 
and consequently the launch region for the wind, grows as 
the SFR increases, then the terminal velocities rise more 
slowly than SFR 5 . Figure 10 illustrates the constant 
surface brightness case previously suggested by Hcckman 
etal. (2000), which implies v te rm oc SFR 1 / 4 . This scaling 
for the size of the launch region fails to reproduce the steep 
slope of the outflow velocity envelope from dwarf starburst 
to LIG luminosities. 

It is not obvious that the hot wind would naturally 
tune the terminal cloud velocities to match the rotation 
speed of the parent galaxy as suggested by Figure 8. This 
correlation implicitly requires p w to be connected to the 
halo mass, presumably via the SFR. Some simple recipes 
are explored here. Using an isothermal halo with ve- 
locity dispersion a to model the halo mass distribution, 
Mir) = 2G~ 1 a 2 r, Equation 5 becomes 



dv c _ -2a 2 
dr 



v r .r 



1 



v r r* 



(9) 



which is easily solved for v(r) when v c << v w . In or- 
der for the clouds to move out a significant distance, the 
solution suggests the characteristic flow velocity v c har = 
\J"irri w v w I (2QNh Rom p [i) must exceed 2a (e.g. Mur- 
ray, Quataert, & Thompson 2004). For purposes of 
illustration, suppose the star formation rate scales as 
SFR cx f g M tot /t dyn cx a 3 , where f g M tot is the gas 
mass of the galaxy. This scaling is plausible since 
the mass of a galaxy increases as a 3 and the dynami- 
cal timescale of a disk, R^/V, is independent of mass 
(Mo, Mao, White 1998). Normalizing to the stellar ve- 
locity dispersions in IRAS23365+3604, IRAS20087-0308, 
and IRAS17208-0014 (Genzel etal. 2001) gives SFR w 
500 M© yr _1 (cr/210 km s -1 ) 3 . Using this relation and in- 
tegrating Equation 9 from low initial velocity at i?o = 
200 pc, clouds accelerate to the velocities shown in the 
top panel of Figure 11 before gravity slows them down. 
Choosing (N H /10 22 cm~ 2 )(i? /200 pc) w 1 again provides 
the best normalization to the data. With this SFR cx a 3 
recipe, the cloud velocities rise steeply with the rotation 
speed because star formation is greatly suppressed in the 
low mass galaxies. Scaling the SFR as a 2 or a produces 
a much more gradual rise in outflow speed with rotation 
speed, and these models are not consistent with the data. 

This exercise demonstrates that ram pressure from a 
supernova-heated wind accelerates cool gas to a maximum 
velocity determined by the temperature of the momentum- 
carrying component of the hot wind. Outflow velocities 
measured for the 2 Jy sample of ULIGs are consistent with 
a maximum wind speed near v 500 km s _1 , but more 
measurements will be required to confirm or refute this 
suggestion. Indeed, in Figure 10, four galaxies from the 
Rupke et al. and Heckman et al. samples present outflows 
significantly faster than this fiducial maximum. Several of 
the 2 Jy ULIGs also show absorption, albeit far from line 
center, at velocities in excess of 500 km s _1 . Failure of 
the outflow measurements to define an asymptotic velocity 

12 For this to be strictly true, the ionization parameter must exceed Uo 
recombination rate and the rate of dust absorption per H atom. ^ > 



might reveal (1) the presence of a hotter wind component 
or (2) an acceleration mechanism other than supernova. 
The awkwardness of the supernova-driven wind model for 
explaining the terminal velocity - rotation speed correla- 
tion motivates some consideration of alternative accelera- 
tion mechanisms. 

5.2. Radiatively- Driven Winds 

Stellar wind velocities exhibit a scaling similar to that 
found here for galactic winds. The terminal wind velocities 
from central stars of planetary nebulae are typically 4.4 
times the escape velocity at the stellar surface (Lamers & 
Cassinelli 1999). Is this remarkable coincidence an indica- 
tion that the underlying wind physics is similar? The cor- 
relation seems natural in stellar winds. The wind is accel- 
erated by radiation pressure on dust grains, and the max- 
imum surface luminosity is limited by the requirement of 
hydrostatic equilibrium. In galaxies, the SFR determines 
the momentum of the hot wind, but the relation between 
galaxy mass and SFR is complicated by variable quanti- 
ties such as the gas fraction and the star formation effi- 
ciency. Radiatively-driven winds on a galactic scale have 
an appealing property, namely that the galactic luminosity 
cannot increase much beyond an effective Eddington-like 
luminosity. 

The radiative force exerted on the cool gas from dust 
absorption usually dominates over that from photoelec- 
tric absorption by gas. 12 For the radiatively-driven wind 
concept to be viable for galaxies, the dust grains need to 
survive in the outflow and pull the gas along. Dust grains 
in the interstellar medium are charged and typically cou- 
ple electrostatically to the gas (Draine 2003). Sputtering 
following the passage of a shock destroys grains, but calcu- 
lations suggest significant amounts of dust survive in radio 
galaxies (De Young 1998; Villar-Martin etal. 2001). Sev- 
eral types of observations strongly suggest that starburst 
winds are indeed dusty. First, the Na I line is a resonance 
line, so the absorbed photons are re-emitted in random di- 
rections. The absence of Na I emission in the ULIG spec- 
tra presented in this paper suggests the re-emitted photons 
are absorbed by dust (and converted to heat) before they 
can escape the galaxy. Second, the thermal emission from 
cold dust in starburst galaxies has a significantly larger 
scale-height than the stars indicating dust is lofted above 
the region where it forms (Engelbract et al. 2004, NGC 55; 
Popescu etal. 2004, NGC 891). Third, extended red emis- 
sion, a broad emission band seen in many dusty astrophys- 
ical objects, has been detected in the halo of the superwind 
galaxy M82 (Perrin, Darbon, Sivan 1995; Gordon, Witt, & 
Fricdmann 1998). Furthermore, recent GALEX imaging of 
M82 show ultraviolet halo emission significantly brighter 
than expected from the nebular continuum, which likely 
indicates the light is scattered by small grains (Chris Mar- 
tin, pvt comm). The dynamical implications of dust in 
these winds needs to be considered. 

Little ultraviolet radiation emerges from ULIGs imply- 
ing these galaxies are optically thick to much of the in- 
trinsic starburst luminosity. The radiative force on the 
gas clouds is F ra d — M c c~ x J n v F v dv, where the absorp- 
tion opacity of the dust and gas mixture, k, drops from 

0.01 (Dopita et al. 2002). This is easily seen by equating the hydrogen 
£LB_ w here 5 is the photon flux. 

kc euv 
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600 cm 2 g" 1 near 1000 A to 100 cm 2 g _1 in the U band 
(Li & Draine 2001). In the optically thick limit, the radia- 
tion force reaches L jc and balances the gravitational force 
on the cool interstellar gas when the luminosity reaches 
.ed = 4G 1 c/ g cr 4 . For a fiducial halo velocity dis- 
persion of 150 km s _1 , the critical luminosity required to 
radiatively accelerate most of the cool interstellar gas is 
then 

L — - 2 - 4x 1012 L °( iso^a) (10) 

The fraction of the galactic mass in the cool wind, f g , will 
be less than 0.1, so the ULIG luminosities appear high 
enough to radiatively accelerate significant amounts of cool 
gas and dust. Assuming a starburst duration of 100 Myr 
(Lcitherer 1999) and a Salpctcr (a = 2.35) initial mass 
function from 0.1 to 100 M Q , the corresponding star for- 
mation rate where radiation pressure becomes important 
for wind acceleration is 

S F K « rt ,~4OOM 0jt -. (jin^-T) 4 5T <"> 

Any additional (besides supernova) acceleration mech- 
anism should predict the highest outflow speeds and/or 
explain why most outflows are tuned to the rotation (or 
escape) speed. The latter is clear. Once the starburst lu- 
minosity exceeds L cr u t ED, the cloud is accelerated from 
rest at Rq to v(r) = 2o (L j L crit . ED - l)ln(r/R ) (e.g. 
Murray, Quartaert, and Thompson 2004). The weak radial 
dependence indicates the terminal cloud velocity depends 
mainly on the ratio L / L C rit,ED- Figure 11 (bottom panel) 
demonstrates that L/L crit .ED ~ 2 provides an excellent 
description of the correlation between outflow velocities 
and rotation speed. 

The outflow speeds exceeding 500 km s" 1 could have 
several explanations. If radiation pressure is the domi- 
nant acceleration mechanism, then starbursts are not the 
source of luminosity because extrapolation of the velocities 
in Figure 11 to 1500 km s _1 implies rotation speeds larger 
than galaxies. AGN would have to supply the radiative ac- 
celeration. Ram-pressure could also accelerate cool gas to 
these high velocities if the clouds are launched from small 
radii associated with AGN rather than starburst radii. In- 
spection of the starburst spectra with outflows faster than 
500 km s- 1 indicates IRAS-FSC05024-1941 (1538 km s" 1 
, Rupke etal. 2002) has cool infrared colors but is op- 
tically identified as a Seyfert 2; IRAS1 11 19+3257 (1410 
km s _1 , Heckman etal. 2000) has warm infrared color. 
It therefore seems quite likely that both these objects have 
active nuclei, and one can speculate that AGN are respon- 
sible for the abnormally high velocities. A much larger 
database will clearly be required to determine the relative 
roles of AGN and starbursts in driving winds. Among the 
the 2 Jy CO-subsample, however, only IRAS11506+1331 
shows Seyfert 2 activity, and it does not present an unusu- 
ally fast wind. 

6. CONCLUSIONS 

Doppler shifts measured from interstellar absorption 
lines in starburst spectra reveal cool gasesous outflows. 
The fraction of starbursts presenting outflows increases 
with starburst luminosity reflecting either the brief duty 



cycle of the ultraluminous phase or the large solid angle 
subtended by these winds. The cool outflow does not en- 
tirely cover the optical nucleus indicating clouds or shells 
are a better description than a continuous fluid. The cov- 
ering fraction approaches unity in the galaxies with the 
strongest winds (i.e. largest equivalent width), and these 
merging systems are usually found to be near perigalac- 
tic passage. More widely separated systems and single- 
nuclei systems generally have lower covering factors and 
line strengths. These results lead to the conclusions that 
(1) much of the mass loading occurs very near perigalactic 
passage and (2) the mass loading is gradually shut down 
as the nuclei finish merging. The mass entrainment rate 
of cold clouds is quite large, possibly comparable to the 
SFR, and is discussed in more detail in Paper II. 

Galaxies with higher SFRs accelerate the absorbing 
clouds to higher velocities. Projection effects introduce a 
large range in measured velocities at a given SFR, but the 
maximum outflow velocity rises roug hly as SFR 1 / 3 . This 
trend challenges the hypothesis that the clouds are accel- 
erated by the hot winds, whose emission-measure weighted 
temperatures appear to vary little with circular velocity. 
The simplest explanation is that the clouds in the ultra- 
luminous starbursts reach terminal velocities comparable 
to the velocity of the hot winds, and that the low Na I 
wind velocities in dwarf starbursts reflect a shortage of 
momentum in dwarf starburst winds. 

However, in order for momentum-driven winds to ex- 
plain the intriguing similarity between the terminal cloud 
velocity and the galactic escape velocity, the rate of mo- 
mentum injection must be tightly correlated with the 
depth of the gravitational potential. Radiativcly-driven 
winds naturally provide this connection via a galac- 
tic equivalent to the stellar Eddington luminosity, while 
rather ad hoc recipes between the SFR and dark halo 
mass must be cooked up for supernova-driven winds. A 
highly simplified description of how dusty clouds would be 
accelerated by the starburst continuum luminosity demon- 
strates that ULIG luminosities are sufficient to accelerate 
significant columns of cool gas to the observed Na I veloc- 
ities. 

These results draw attention to the importance of feed- 
back in very massive galaxies. Preliminary estimates of the 
mass loss rates suggest a significant fraction of the ISM 
participates in the cool outflows. This result is surpris- 
ing in the context of energy-driven winds because a large 
fraction of the supernova energy is expected to be radi- 
ated away in the dense molecular ISM of an ULIG. Much 
as supernova feedback is most effective in dwarf galaxies 
where hot gas can escape from the halo, radiative feed- 
back is naturally strongest in the densest, most luminous 
galaxies. Hence, radiative feedback may help explain the 
differences between the shape of the halo mass function 
and the galaxy luminosity function at the high-mass end. 
The dynamical properties of the stars in ULIGs do appear 
to be consistent with those of field ellipticals (Genzel et al. 
2001), so it is also appealing to associate the winds seen in 
ULIGs with the event that shuts the star formation down 
quickly enough to create high ratios of alpha elements to 
iron in elliptical galaxies. 

Filling in Figure 10 and Figure 11 with more observa- 
tions will help clarify some of these isssues. Also, since 
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radiatively-driven winds do not require a hot phase, X- 
ray observations of these ULIGs may distinguish acceler- 
ation mechanisms. Spitzer and GALEX observations will 
be instrumental for constraining the amount of extrapla- 
nar dust in starbursts and its composition. A few of the 
most luminous objects present velocities too high to ex- 
plain with starburst winds, so it will be possible to ex- 
plore the nature of outflows in warm ULIGs during the 
SB/AGN transition. More detailed dynamical modeling 
of the survivability and fate of the cool clouds as well as 
a more realistic treatment of radiation pressure will allow 
a more detailed comparison to the data (e.g. line widths, 
spatial gradients, etc.) and thereby a more stringent test 
of the ideas presented here. 



The observations presented here would not have been 
possible without the R band images of the galaxies pro- 
vided by Tom Murphy. The author thanks Lee Armus 
for help with some of the spectral fitting and enthusiasm 
for starting this project. Discussions with Norm Murray, 
Todd Thompson, and Eliot Quataert regarding the theo- 
retical interpretation of the data presented here are great- 
fully acknowledged. This paper was completed at the As- 
pen Center for Physics where discussions with Karl Gor- 
don and Timothy Hcckman improved the content of this 
paper. Funding for this project was provided by the David 
and Lucile Packard Foundation and the Alfred P. Sloan 
Foundation. Finally, the author thanks those of Hawai- 
ian ancestry on whose sacred mountain these data were 
obtained. 
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Table 1 
Sample of Galaxies 



Object 


z 




Lfir 




SFR 




(CO) 


(Lo) 


(Lo) 


( M yr" 1 ) 


IRAS00153+5454 


0, 


1120 


11 


,94 


12 


.27 


320 


IRAS00188-0856 


0. 


1285 


12 


.18 


12 


.41 


450 


IRAS00262+4251 


0. 


.09724 


11 


.90 


12 


.16 


250 


IRAS03158+4227 


0. 


1344 


12 


.39 


12 


.64 


750 


IRAS03521+0028 


0. 


1519 


12 


.33 


12 


.57 


640 


IRAS08030+5243 


0. 


.08350 


11 


.82 


12 


.05 


190 


IRAS 10494+4424 


0. 


09231 


11 


.99 


12 


.23 


290 


IRAS10565+2448 


0. 


04311 


11 


.81 


12 


.05 


190 


IRAS11506+1331 


0. 


.1273 


12 


.11 


12 


.35 


390 


IRAS15245+1019 


0, 


.0757 


11 


.81 


12 


.10 


220 


IRAS16090-0139 


0. 


1336 


12 


.34 


12 


.57 


630 


IRAS 16487+5447 


0, 


1036 


11 


.98 


12 


.29 


330 


IRAS17208-0014 


0. 


.04282 


12 


.23 


12 


.45 


490 


IRAS18368+3549 


0. 


.1162 


12 


.03 


12 


.27 


320 


IRAS19297-0406 


0. 


.08573 


12 


.21 


12 


.44 


470 


IRAS19458+0944 


0. 


1000 


12 


.15 


12 


.39 


420 


IRAS20087-0308 


0. 


1057 


12 


.23 


12 


.47 


510 


IRAS23365+3604 


0. 


.06448 


11 


.96 


12 


.20 


280 



Table Notes: (1) Object name. (2) Redshifts from Solomon et al. (1997) except IRAS 00153+5454, IRAS 15245+1019, 
and IRAS 16487+5447 which are from Dr. Aaron Evans (pvt. comm.). Velocities have errors of ±cz w 20 km s^ 1 . (3) 
The far-infrared luminosity computed from Lfir = 3.86 x 10 5 L d| / (2.58F !y (6O?im) + f v (100um)), where the flux density 
is in Janskys and the luminosity distance (h=0.7, Sl = 0-3 and J7a = 0.7) is in Mpc. (4) The estimated 8 - 1000 um 
luminosity Lir (h=0.7, Qq = 0.3 and J7a = 0.7). The IRAS 12um fluxes for the galaxies in this sample are upper limits, 
so the median flux ratios of bright ULIGs with IRAS detections in all four bands were used to estimate Lir as described 
in Murphy etal. (1996). (5) The star formation rate estimated from SFR = Lir/5.8 x 10 9 L Q (Kennicutt 1989), where 
the luminosity is the estimated bolometric luminosity Lir. This relation uses the continuous star formation model of 
Lcitherer etal. (1999) and a Salpeter initial mass function from 0.1 to 100 M Q . 
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Table 2 

Properties of Na I Absorption 



Object 


IF (NaD) 


/. 


Umax 


Cf 


Al 




(A) 


(A) 


( km s" 1 ) 






(Vo 1 k P c ) 


IRAS00153+5454 


4.1 ±0.5 


0.15 


-715 





.14 


3.8 (773) 


IRAS00188-0856 


7.5 ±0.6 


0.09 


-751 





.19 


4.6 (179) 


IRAS00262+4251 


1.0 ±0.3 


0.33 


-580 





.11 


6.9 (378) 


IRAS03158+4227 


5.1 ±0.4 


0.052 


-1047 





.25 


6.2 (276) 


IRAS03521+0028 


5.1 ±0.4 


0.1 


-350 





.63 


5.0 (179) 


IRAS08030+5243 


7.4 ± 1.0 


0.063 


-311 





.26 


> 8.2 (572) 


IRAS10494+4424 


5.6 ±0.1 


0.024 


-560 





.15 


> 12.0 (770) 


IRAS10565+2448 


10.4 ±0.1 


0.03 


-620 





.42 


> 11.4 (1374) 


IRAS11506+1331 


3.88 ±0.2 


0.10 


-435 





.28 


4.53 (179) 


IRAS15245+1019 


8.5 ±0.1 


0.05 


-559 





.67 


4.2 (1274) 


IRAS16090-0139 


3.32 ±0.06 


0.083 


-543 





.11 


7.6 (372) 


IRAS 16487+5447 


< 0.5 


> 0.46 








2.71 (577) 


IRAS17208-0014 


9.4 ±0.1 


0.06 


-691 





.52 


> 7.3 (876) 


IRAS18368+3549 


9.0 ±0.3 


0.04 


-768 





.25 


10.3 (478) 


IRAS19297-0406 


3.7 ±0.2 


0.12 


-697 





.22 


4.8 (279) 


IRAS19458+0944 


4.0 










IRAS20087-0308 


9.0 ± 1.5 


0.074 


-1082 





.26 


12.2 (673) 


IRAS23365+3604 


4.5 ±0.4 


0.12 


-651 





.16 


5.9 (477) 



Tabic Notes - 

(1) Object name. 

(2) Total measured NaD absorption equivalent width (observed frame). 

(3) Fraction of total NaD equivalent width from a stellar component predicted from measured Mgb equivalent width and 
the relation VF(Na I) = l/3W(Mg I) described in the text. 

(4) Terminal velocity measured from the intersection of the bluest part of the absorption line with the continuum. 

(5) Fraction of continuum source covered by absorbing clouds where C/ = 1 — /(A5890). 



15 



Table 3 

Fitted Na I Absorption (Optically Thick) 





W Don 


V B 




An 








(k) 


( \tyt\ c ^ 


( Irm c ^ 


( \r m o ^ ~\ 
I JSTll O J 




Cm 12 rm- 2 \ 


IJxAoUU 100+0404 


O. / 


— zyz zt 10 


a -U /in 

40Z zt 4U 


KKA -U 01 O 

oo4 zt ziy 


^ l.UD ± U.UD 


436_ 392 


IRAS00188-0856 


4.42 


-328 ±30 


520 ±20 


260 ± 11 


1.07 ±0.09 


310t 2 f 


IRAS00262+4251 


0.94 


-293 ±20 


296 ± 20 


500 ± 150 


1.7±0.4 


1 ro+0.94 

l.OO_ 25 


IRAS03158+4227 


4.43 


-499 ±25 


405 ± 35 


312 ±90 


1.21 ±0.15 


251* 2 


IRAS03521+0028 


5.0 


-77 ± 12 


175 ±31 


636 ± 360 


1.42 ±0.46 


l.67%\ 


IRAS08030+5243 


4.02 


+267 ±2 


158 ±4 


350 ± 30 


1.16 ±0.05 


29-8t^ 5 


IRAS 10494+4424 


2.85 


-270 ±23 


315 ±35 


237 ±7 


1.07±0.18 


174ii? f o 


IRAS10565+2448a 


7.80 


-307+ 16 


307 ± 16 


278 ± 107 


< i.io 


> 160 


IRAS10565+2448b 


5.1 


-87 ±2 


188 ± 10 


11 


1.33 ±0.18 


19.6 20 4 6 3 


IRAS11506+1331 


3.3 


-96 ±6 


242 ± 10 


222 ± 150 


1.10 ±0.07 


59.8t^ f 4 


IRAS15245+1019 


8.9 


-240 ± 4 


301 ±6 


347 ± 108 


1.24 ±0.03 


43.ltf. 2 


IRAS16090-0139 


1.31 


-282 ±23 


296 ± 25 


354 ±22 


1.96 ±0.28 


z - i -0.2 


IRAS16487+5447 












IRAS17208-0014 


6.91 


-409 ± 3 


292 ±4 


254 ±6 


1.06 ±0.02 


883±£ f 2 


IRAS18368+3549a 


4.64 


-420 ± 35 


371 ± 49 


364 ± 134 


1.06 ±0.28 


5671$ 


IRAS18368+3549b 


4.89 


-143 ±4 


213 ± 12 


11 


1.24 ±0.24 


22.6+^ 


IRAS19297-0406 


3.84 


-114 ± 16 


377 ± 26 


313 ±60 


1.46 ±0.18 


11+5.4' 
iJ --1.9 


IRAS19458+0944 












IRAS20087-0308 


4.3 


-420 ± 50 


628 ±57 


478 ± 113 


1.87 ±0.36 


O or+2.95 
y-^O-0.44 


IRAS23365+3604 


1.80 


-389 ± 72 


308 ± 85 


610 ± 195 


1.26 ±0.79 


o 97+ inf 
°- z ' -6.04 



Table Notes - 

(1) Object name. Spectral model is the sum of a Na I doublet at the systemic velocity and a doublet, the dynamic 
component, for which the velocity is fitted. The doublet ratio, D = 1T(A5890)/1T(A5896), of the systemic component was 
held at D = 1.06 - equivalent to an A5890 optical depth of 142. The doublet ratio of the dynamic component was fitted. 

(2) Equivalent width of fitted Doppler-shifted component. 

(3) Fitted velocity of Doppler shifted interstellar component. 

(4) FWHM of fitted Doppler shifted interstellar component. 

(5) FWHM of the component at the systemic velocity. Parentheses indicate an assumed, rather than fitted, value. 

(6) Fitted doublet ratio, D = 1T(5895)/1T(5889) of Doppler shifted interstellar component. 

(7) Constraints on column density from curve of growth analysis and doublet ratio from column 3. 
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Table 4 

Fitted Na I Absorption (Optically Thin) 



Object 


W Dop 


V B 


Av B 








N NaI 




(A) 


( km s -1 ) 


( km i 


) 


( km t 


r 1 ) 


(10 12 cm" 2 ) 


IRAS00153+5454 


3 


.2 


-289 + 42 


499 


± 


48 


293 


± 


42 


9.94 ±0.67 


IRAS00188-0856 


4. 


.3 


-340 ± 28 


567 


± 


36 


332 


± 


15 


13.1 ±0.6 


IRAS00262+4251 


0, 


.5 


-406 ± 20 


232 


± 


21 


493 


± 


46 


1.5 ±0.6 


IRAS03158+4227 


5, 


.3 


-469+ 11 


466 


± 


29 


340 


± 


83 


15.67±0.33 


IRAS03521+0028 


3. 


.8 


-94+15 


152 


± 


46 


530 


± 


200 


13.1 ±4.3 


IRAS08030+5243 


4. 


.2 


+268 + 5 


166 


± 


16 


536 


± 


397 


13.1 ±1.1 


IRAS 10494+4424 


1 


.2 


-380 + 57 


356 


± 


80 


318 


± 


15 


3.84 ±0.65 


IRAS10565+2448a 


7. 


.7 


-252 + 22 


469 


± 


24 


283 


± 


29 


25.14±0.11 


IRAS10565+2448b 


3 


.9 


-88 + 3 


156 


± 


5 


11 






12.94 ±0.04 


IRAS11506+1331 


3 


.4 


-58 + 2 


338 


± 


27 


305 


± 


90 


10.37 ±0.02 


IRAS15245+1019 


7. 


,4 


-243 ± 6 


296 


± 


8 


307 


± 


20 


23.31 ±0.04 


IRAS16090-0139 


0. 


.7 


-325 + 40 


281 


± 


37 


479 


± 


29 


2.16 ±0.36 


IRAS 16487+5447 






















IRAS 17208-0014 


7. 


4 


-359+ 15 


462 


± 


38 


253 


± 


13 


21.18 ±0.08 


IRAS18368+3549a 


5, 


.4 


-341 + 56 


469 


± 


66 


354 


± 


37 


16.59 ±0.17 


IRAS18368+3549b 


3 


.0 


-140 + 7 


168 


± 


19 


11 






9.46 ±0.10 


IRAS 19297-0406 


3. 


.6 


-413+13 


370 


± 


22 


286 


± 


36 


11.3±0.6 


IRAS19458+0944 






















IRAS20087-0308 


4. 


.6 


-421 + 50 


631 


± 


75 


525 


± 


703 


14.3 ±0.8 


IRAS23365+3604 


1 


,3 


-384 + 96 


302 


± 


22 


660 


± 


179 


4.3 ± 1.5 



Table Notes - 

(1) Object name. Spectral model is the sum of a Na I doublet at the systemic velocity and a doublet, the dyn 
component, for which the velocity is fitted. The doublet ratio was fixed at the optically thin limit. 

(2) Equivalent width of fitted Doppler-shifted component. 

(3) Fitted velocity of Doppler shifted interstellar component. 

(4) FWHM of fitted Doppler shifted interstellar component. 

(5) FWHM of the component at the systemic velocity. Parentheses indicate an assumed, rather than fitted, value. 

(6) Column density of Na I. 



18 




Fig. 1. — (b) The ESI spectra around the Na I A5890, 96 doublet. The emission line is He I A5876. No Na I absorption was detected in 
IRAS 16487+5447 and IRAS19458+0944. The dotted line shows a doublet at the systemic velocity, which was determined from observations of 
the molecular gas. The doublet ratio of the systemic component is fixed at EW (5890) / EW (5896) = 1.06 here. The Doppler shift and doublet 
ratio of the dynamic component, shown by the thin solid line, were fitted. The heavy solid line shows the sum of the dynamic component, 
systemic component, and the He I emission line. All but one galaxy presents a highly blueshifted dynamic component. 
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Fig. 2. — Measured Na D absorption equivalent width versus measured Mg I absorption equivalent width in ULIGs. The solid line is a fit 
to stellar spectra (see §3.1 of text), and the dotted line illustrates typical values in elliptical galaxies (Bica ct al. 1991). 
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-1000-500 500 1000 
V (km/s) 

Fig. 3. — The K I 7664.9 profile for the five galaxies with the strongest interstellar absorption. The weaker member of the doublet is at 
+1327 km s" 1 (not shown). The Na I doublet, arbitrarily scaled, is overlaid. The profiles deviate, as they must, redward of the systemic 
velocity where the Na I 5890 appears. The spikes in the K I spectra are residuals from the subtraction of bright night sky lines and should 
be ignored. The K I 7664.9 absorption profile is generally similar in shape to the Na I 5889 profile but has poor signal-to-noise ratio. 
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Fig. 4. — Examples of fitted Na I absorption lines in the optically thin (top panel) and thick (bottom panel) limits. The blue line shows the 
fitted dynamic component, and the doublet at the systemic velocity is drawn in black. The green line shows the sum of the model components. 
The optically thick models describe the local minima better than the optically thin models. 
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Fig. 5. — (Left ) The Na I equivalent width shows no correlation with the velocity width of the line. (Right) Cloud covering factor versus 
Na I equivalent width in the dynamic component. The sightlines with the highest covering factors present the strongest absorption lines. 
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Fig. 6. — Outflow velocities at line center vs galactic star formation rate. The new data, 1 Jy ULIGs (Rupke et al. 2002), LIGs (Heckman 
etal. 2000), and dwarf starbursts (Schwartz & Martin 2004) are denoted by squares, hexagons, circles, and triangles, respectively. These 
star burst systems represent the highest surface brightness objects in their luminosity class and are thought to define the maximum outflow 
velocities reached. The dashed line represents the fitted upper envelope when a simple model for projection effects is applied (see text in §4.2 
for detailsR). 



24 



> 

> 

cti 



> 



3 



2 - 








1 2 

log SFR (Mo/yr) 



3 



8 - 



6 



4 - 








v 



avg v 



Fig. 7. — (top) Distribution of outflow velocities normalized by the mean at a given star formation rate, (bottom) Histogram of top panel 
for all values of the SFR. The expected number of galaxies with outflow velocity between v n and v m is N(v n ,v m ) = N(v ma x , 0) 2")<^> 
illustrated in the bottom panel by a dashed line. The four objects with v >> v max present a significant departure from the model. 
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Fig. 
speed. 
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i. — (top) Outflow velocity at Na I line center versus galactic rotation speed, (bottom) Terminal Na I velocity versus galactic rotation 



Fig. 9. — Wind strength at progressive stages of the merger. The absorption lines are strongest when the merging galaxies are near 
perigalacticon. They become weaker around apogalacticon and again following the merger of the nuclei. 
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Fig. 10. — Terminal velocity of interstellar gas clouds accelerated by a hot wind. For consistency with X-ray measurements, the hot wind 
velocity is taken as 500 km s _1 , independent of galaxy mass; and the mass flux in the wind is assumed equal to the SFR. Solid lines illustrate 
a model where the clouds are launched at a radius of 200 pc. Dotted lines show the effect of increasing the size of the launch region with the 
square root of the SFR as suggested by Heckman et al. (2000). The results illustrate the critical SFR below which the hot wind momentum 
is insufficient to accelerate the clouds to the velocity of the hot wind. 
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Fig. 11. — Models for the terminal outflow velocity with respect to the galactic rotation speed, (top) The SFR increases as the cube, square, 
and first power of the halo velocity dispersion; and momentum from the associated supernovae accelerates cool clouds to the velocity of the hot 
wind. The models were normalized to the ULIG data such that SFR = 500 M Q yr -1 (cr/210 km/s) 3 , SFR = 500 M yr- 1 (cr/225 km/s) 2 , 
and SFR = 500 Mq yr _1 (a/225 km/s). (bottom) The luminosity is two or three times the Eddington luminosity where L^rji = 4G~ 1 <t 4 c/ c , 
and the clouds, assumed to be optically thick at a radius ten times the launch radius, are radiatively accelerated. The radiative acceleration 
model describes the data as well, and perhaps even better, than than the traditional supernova-driven wind model. 
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